ARTERIAL SPIN LABELING USING TIME VARYING GRADIENTS 

This application is a continuation of prior application serial number 10/121,963, 
filed April 12, 2002, entitled ARTERIAL SPIN LABELING USING TIME VARYING 
5 GRADIENTS, which is incorporated herein by reference. 

Field Of The Invention 

The present application relates to a technique for imaging fluid flow 
noninvasively using MRI with radio frequency arterial spin labeling and, more 
10 particularly, to a labeling technique that permits extension of this technique to a multi- 
slice or volume examination. 

Background Of The Invention 

Magnetic Resonance Imaging (MRI) is an imaging technique used primarily in 
15 medical settings to produce high quality images of the inside of the human body. MRI is 
based on principles of nuclear magnetic resonance (NMR), a spectroscopic technique 
used to obtain microscopic chemical and physical information about molecules. 

MR (magnetic resonance) techniques may be used to study fluid flow, such as, 
for example, blood flow and blood perfusion in tissue. One of many possible 
20 applications is to study blood perfusion in the human brain. 

Typically, during MR studies, images are taken of a volume of interest or of 
"slices" in a region of interest. A slice may be defined as a relatively thin region, which 
may be imaged in a single image. Each slice is said to be composed of several volume 
elements called voxels. Traditionally, volume of a voxel may be approximately 3mm A 3, 
25 however, it may vary depending on the application. The magnetic resonance image is 
said to be composed of several picture elements called pixels. The intensity of a pixel is 
proportional to the NMR signal intensity of the contents of the corresponding volume 
element, or voxel, of the object being measured. Present-day MR studies typically call 
for more than one slice image to be taken, or one or more images of the volume of the 
30 region of interest to be taken. 

Magnetic resonance imaging is based on the absorption and emission of energy in 
the radio frequency range of the electromagnetic spectrum. The human body, for 
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example, is composed in large part of water, which has many hydrogen atoms. 
Hydrogen nuclei have an NMR signal in the presence of the magnetic field, and after 
being perturbed by RF (radiofrequency) irradiation. Protons in the hydrogen atoms (and 
in other atoms) possess a property called spin, which can be thought of as a small 
5 magnetic field and which will cause the nucleus to produce an NMR signal. Water 
consists in large part of hydrogen atoms. There are known apparatuses for detecting 
water in tissue using the NMR signal of hydrogen nuclei. Using these basic principles, 
fluid content may be measured in a variety of substances or tissues. 

In order to detect fluid flow or perfusion in a particular region of interest, fluid 
10 flowing into that region may be "labeled" by reversing, or perturbing, the spins of the 
protons of the fluid in some region that is "upstream" from the region of interest, and 
then detecting the labeled fluid when it flows through or is perfused in the region of 
interest. 

Such studies have previously been performed by applying a constant magnetic 
15 field to the fluid and applying RF irradiation to label the spins. In order to spatially 
isolate labeling areas, a magnetic field gradient may be used. Spatially isolating a 
particular area means applying, in the particular area, a magnetic field different in 
strength from the magnetic field in other areas. By using the magnetic field gradient to 
spatially isolate particular areas, only fluid in a predetermined region may be labeled. 
20 However, magnetization transfer effects and other unrelated errors may interfere with the 
procedure by causing more than just the "upstream" fluid to be labeled. In order to 
account for such effects, a control procedure may be used, where a magnetic field 
gradient and RF irradiation are applied so as not to label the spins, but to mimic the 
unrelated effects. 

25 One of the concerns in performing such a study is possible power deposition 

resulting from application of RF irradiation. Typically, power deposition is higher where 
the present magnetic field and the intensity of RF irradiation are higher. There are 
boundaries on the power that may be deposited in the human body without harming the 
body. However, stronger RF irradiation results in higher sensitivity that in turn produces 

30 better results. Therefore, there is a need for a higher-sensitivity method for performing 
fluid flow studies, especially multi-slice or volume studies, in the presence of strong 
magnetic fields, while taking into account maximum power deposition considerations. 
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Summarv Of The Invention 

The present invention relates to a method and apparatus for conducting a 
magnetic resonance fluid flow study by conducting a labeling procedure and a control 
5 procedure and combining datasets from those procedures to create a dataset for the fluid 
flow study. An amplitude modulated magnetic field gradient and amplitude modulated 
RF irradiation may be applied during a labeling procedure. Likewise, amplitude 
modulated RF irradiation may be applied during control procedure. An amplitude 
modulated magnetic field gradient may also be applied during the control procedure. 

10 An envelope for modulating the label RF irradiation may be an absolute value of 

an envelope for modulating the control RF irradiation, such that an average amplitude for 
the label RF irradiation is not zero. The average amplitude for the label RF irradiation 
may be positive or negative, depending on the requirements of a particular application. 
An average amplitude for control RF irradiation may be zero or near zero. 

15 An amplitude modulated magnetic field gradient may also be applied during the 

control procedure. An envelope for modulating the control magnetic field gradient may 
be similar to the envelope for modulating control RF irradiation. 

The control RF irradiation may further be frequency modulated. An envelope for 
frequency modulation of the control RF irradiation may be similar to the envelope for 

20 amplitude modulation of the control RF irradiation. In one embodiment of the invention, 
the envelope for amplitude modulation of the control RF irradiation may be a modified 
square wave. The envelope for the label RF irradiation may be the absolute value of the 
modified square wave. 

An image of fluid flow may be generated by subtracting the label dataset from 

25 the control dataset. In an alternative embodiment of the invention, the image of the fluid 
flow may further be modified by subtracting artifacts of systematic errors from the 
image. 

Another aspect of the invention is a method for compensating for magnetization 
transfer effects by applying an amplitude modulated magnetic field gradient and 
30 amplitude modulated RF irradiation during a label procedure and applying an amplitude 
modulated control magnetic field gradient and amplitude modulated control RF 
irradiation during a control procedure. Such compensation for the magnetization transfer 



effects may be performed in a context of an MR study, such as, for example, an MR 
imaging study. Furthermore, RF irradiation may be frequency modulated during the 
label and/or control experiments. An envelope for label RF irradiation for frequency and 
amplitude modulation may be similar and may be similar to an envelope for amplitude 
5 modulation for label magnetic field gradient. An envelope for control RF irradiation for 
frequency and amplitude modulation may be similar and may be similar to an envelope 
for amplitude modulation of the control magnetic field gradient. The envelope for label 
RF irradiation modulation may be an absolute value of the envelope for the control RF 
irradiation. The envelope for the control RF irradiation may be a modified square wave. 

10 In an alternative embodiment, the envelope for the control RF irradiation may be a 
function other than the square wave, as long as an average control RF irradiation 
amplitude is at or near zero. 

In yet another aspect of the invention, a magnetic resonance apparatus may be 
constructed and arranged to perform fluid flow imaging. Such apparatus may apply 

15 amplitude modulated label RF irradiation and a magnetic field gradient during the 

labeling procedure. Furthermore, such apparatus may apply amplitude modulated label 
RF irradiation and a magnetic field gradient during the control procedure. The envelope 
for applied control magnetic field gradient may be similar or identical to the envelope for 
applied control RF irradiation. The envelope for applied label RF irradiation and the 

20 magnetic field gradient may be the absolute value of the envelope for the control RF 
irradiation. 

In yet another embodiment of the invention, a computer readable medium may be 
provided. Encoded on the computer readable medium may be instructions that, when 
executed, direct a method for generating an image of fluid flow obtained from a 

25 magnetic resonance system. Such a method may comprise combining a control image 
received from the control procedure where the amplitude modulated RF irradiation and 
the magnetic field gradient were applied with a label image received from the labeling 
procedure where the amplitude modulated RF irradiation and the magnetic field gradient 
were applied. The control and label images may be combined to create the image of 

30 fluid flow by subtracting the label image from the control image. 

In another embodiment of the invention, the image of the fluid flow may be 
further modified by subtracting artifacts of residual or system errors. A dataset for the 



residual or system errors may be obtained by performing labeling and control procedures 
in a region different from the region being imaged such as upstream or downstream from 
the region being imaged. The fluid flow image may be further modified as appropriate 
for a particular application. 

5 

Brief Description Of The Drawings 

In the drawings, which are incorporated herein by reference, and in which like 
elements have been given like reference characters, 

Fig. la schematically illustrates a conventional single-slice MR blood flow 
10 mapping study; 

Fig. lb schematically illustrates how a conventional magnetization transfer 
compensation scheme fails to compensate for magnetization transfer effects when used 
in a multislice MR study; 

Fig. 2 schematically illustrates blood flow after spin inversion and a Bloch 
15 equation for computing spin relaxation effects; 

Fig. 3 illustrates a pulse sequence timing diagram for spin labeling and image 
acquisition; 

Fig. 4 is a flowchart illustrating one embodiment of the invention; 
Fig. 5 is a diagram of an amplitude modulation envelope for labeling and control 
20 RF irradiation and for control and labeling magnetic field gradient; 

Fig. 6 is a diagram of a frequency modulation envelope for labeling and control 
RF irradiation; 

Fig. 7 illustrates a direction of magnetic fields with simultaneous RF and 
magnetic field gradient amplitude modulation; 
25 Fig. 8 is a diagram of simulated residual subtraction error; 

Fig. 9 is a diagram of a simulated efficiency of the first embodiment of the 
invention; 

Fig. 10 is an illustration of in- vivo perfusion imaging achieved with the first 
embodiment of the invention; and 
30 Fig. 1 1 is an illustration of comparative results of labeling inferior to an imaged 

region and superior to the imaged region. 



-6- 



Detailed Description 

One embodiment of the present invention is described below in relation to Figs. 
2-11. Although the embodiment of Figs. 2-1 1 refers to measuring perfusion of blood in 
the brain, the present invention is not so limited. The invention may also be applied to 
5 any MR study of fluid flow or perfusion, whether in live tissue, or under other 

conditions. Furthermore, the present invention may be used for angiographic studies and 
other medical and scientific applications, as deemed appropriate by one skilled in the art. 

The embodiment of the invention illustrated in Figs. 2-1 1 involves proton 
magnetic resonance imaging of perfusion using water as a freely diffusable tracer. 

10 Application of the invention to the measurement of cerebral blood flow (CBF) is detailed 
in relation to figures 10 and 1 1 . 

Although terms "flow" and "perfusion" may sometimes be used interchangeably, 
perfusion as used herein refers to a diffusable exchange between a fluid and a substance, 
such as, for example, human tissue. The term "flow" as used herein, generally refers to 

15 flow of liquid in vessels, such as, for example, flow of blood in arteries. 

Generally, the method of measuring blood flow or perfusion involves labeling 
proton spins of inflowing water in the arterial blood using magnetic resonance. The 
blood or water containing hydrogen atoms with spins exchange with bulk water in the 
brain eventually reaching a state which exhibits a decreased intensity in a brain image, 

20 which may be detected using MR imaging techniques. This method of measuring fluid 
perfusion is generally referred to as Arterial Spin Labeling (ASL). Furthermore, such a 
method may be performed by applying RF irradiation to a certain volume of blood to 
label it - such a method is more specifically referred to as Continuous Spin Labeling 
(CASL). One aspect of the present invention relates to compensating for magnetization 

25 transfer effects in a multi-slice or volume study using CASL. 

In some arterial spin labeling methods, the labeling is achieved with a brief pulse 
of RF such that the position of the flowing spins does not change significantly compared 
to the smallest dimension of the region where the labeling has its greatest effect. These 
methods are hereafter referred to as "pulsed" ASL. The duration of such pulsed labeling 

30 may depend upon the application but in the application of measuring perfusion to the 
brain, it is typically much less than 100ms. Other methods for ASL apply RF over a 
period of time such that spins can readily flow in and out of the region where the labeling 



has its greatest effect. Often the motion of the fluid is included in the design and 
optimization of the labeling method. Such ASL methods are hereafter referred to as 
continuous ASL. The term "continuous" does not preclude modulation of amplitude or 
another property of the RF or other fields on a time scale faster than the total duration of 
the RF waveform which helps to achieve the labeling of the flowing fluid. Nor does the 
term "continuous" preclude a long pause between the last time at which the RF 
amplitude is greater than zero and the measurement of signal. 

For simplicity, in the foregoing description, labeling procedures are described as 
taking place before control procedures. However, the invention is not so limited. 
Multiple control and labeling procedures may be performed. Such multiple procedures 
may be used to generate averaged images, as described below in connection with Figs. 
10 and 1 1 . Labeling and control may be arbitrarily interleaved, as determined by one 
skilled in the art. 

A single-slice ASL imaging technique is now described in connection with 
Figures la and lb. In such a technique, RF irradiation is applied to volume 12 in order 
to label the spins. The direction of blood flow 10 carries the labeled spins to image slice 
14a in the region of interest. In order to compensate for magnetization transfer effects, 
control RF irradiation may be applied to volume 1 6, outside of the region of interest. 
The spins are not labeled in the control procedure, because volume 16 is "downstream" 
from blood flow. In an alternative embodiment, volume 16 may have no blood flow. 
The magnetization transfer effects on slice 14a will be the same in the labeling and 
control procedures only if regions 12 and 16 are equidistant from slice 14a. Both control 
and labeling irradiation may be constant and applied in the presence of the constant 
magnetic field gradient. 

In order to produce an image of the blood perfusion in slice 14a, images of 
control and labeling may be taken and combined to produce the perfusion image. The 
perfusion image may be produced, for example, by subtracting the label image from the 
control image. 

Such a single slice ASL technique may not be applicable to multi-slice studies, as 
illustrated in connection with Figure lb, where imaging of multiple slices 14x is shown 
schematically. In such a study, distance between the labeling region 12 and a slice and 
distance between the control region 16 and the slice will not be equal for all the slices - 



such as, for example, slices 14a and 14c. Therefore, the magnetization transfer effects 
will not be properly accounted for in the control for those slices, and the images received 
from combining of control and labeling images will not be accurate. 

Previously, a different technique was used to perform multi-slice study. In this 
conventional technique, a constant magnetic field gradient and a constant RF irradiation 
may be applied in labeling. In control, however, amplitude modulated RF irradiation 
may be applied in the presence of the constant magnetic field gradient. This amplitude- 
modulated RF irradiation is applied to the same region as labeling irradiation, but, 
because of the amplitude modulation, it generally results in not labeling the spins. The 
sensitivity of this technique in practice may depend on the strength of the RF irradiation 
overall. 

One embodiment of the invention will now be described in connection with 
Figures 2-11. Fig. 2 schematically illustrates blood flow after spin inversion. Spin 
inversion is performed by RF irradiation in inversion plane 12, after which blood flows 
through the arteries 22 proximal to tissue 14 in the region of interest; after which blood is 
perfused in tissue 14. 

Spin inversion is achieved by applying magnetic fields in the region where 
nuclear spins are located. The magnetization, or net magnetic properties of a region 
containing nuclear spins, changes direction and amplitude depending upon the magnetic 
fields in that region and times Tl and T2, material specific properties of the substance 
within the region. The change of the magnetization vector with time is determined by the 
Bloch equations, which are familiar to those skilled in the art. 

These constants - Tl (spin lattice relaxation time) and T2 (spin relaxation time) - 
may be used, along with other known values, such as the speed of the blood flow, to 
calculate the magnetization effects on the blood at a particular time t after labeling. 
Equations describing the behavior of the magnetization vector under various conditions 
are called Bloch equations. 

Equation 24 illustrates behavior of the magnetization vector for the blood flow 
path illustrated in Fig. 2: 
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Control RF irradiation may be performed in the same plane 12, however, it may 
be performed such that the spins are not inverted, while the magnetization transfer effects 
and unrelated errors are the same as in the labeling procedure. 

If imaging is performed at time t illustrated in Fig. 2, the signal from arteries 22 
may overpower the perfusion signal in tissue 14. Furthermore, not all areas of tissue 14 
are covered by labeled blood. Therefore, it may be beneficial in some applications to 
delay imaging after labeling for some period of time. 

Fig. 3 illustrates a possible sequence timing diagram for spin labeling and image 
acquisition. During labeling period 32, labeling RF irradiation is applied to a particular 
region, after which blood flows into the region of interest during delay 34, and, finally, 
imaging is performed during period 36. During period 36 one or more slices may be 
imaged. In an alternative embodiment, volume imaging may be performed and one or 
more volume images may be taken. "Imaging," as used herein, generally refers to taking 
data sets from the region of interest and then manipulating them as appropriate, including 
producing a visual image, if desired. The same or a similar timing diagram may be used 
for control imaging. 

In one embodiment of the invention, labeling RF irradiation may be amplitude 
and frequency modulated. Furthermore, the labeling magnetic field gradient may also be 
amplitude modulated. In order to control for magnetization transfer effects properly, the 
control RF may also be amplitude and frequency modulated and control magnetic field 
gradient may also be amplitude modulated. For both labeling and control, amplitude 
modulation of the RF irradiation, frequency modulation of the labeling RF irradiation, 
and amplitude modulation of the magnetic field gradient may all be performed 
simultaneously. Each modulation may be performed according to a particular 
modulation envelope. The term "envelope" is used herein to refer to a function or a 
curve description according to which a particular parameter may be modulated. The 
envelope may be a general description, and the real modulation envelope may differ from 
application to application, as deemed appropriate by one skilled in the art. Furthermore, 
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the particular scale of modulation may depend on the kind of modulation and may differ 

in each particular implementation. 

Fig. 4 is a flowchart illustrating the steps carried out in accordance with one 

embodiment of the invention. A subject of the MR study and the apparatus are prepared 
5 in step 402, after which an amplitude-modulated magnetic field gradient may be applied. 

Labeling amplitude modulated RF irradiation may be applied to perform spin inversion. 

Amplitude modulation may be performed according to a labeling modulation envelope. 

Labeling RF may also be frequency-modulated according to a labeling envelope. 

After labeling RF irradiation is applied and a passage of a predetermined delay 
10 period (Fig. 3), a first set of MR data may be acquired in the region of interest in step 

406. Such a dataset may contain data for MR imaging or other relevant data, as 

determined by one skilled in the art. 

Residual effects of spin labeling may be erased as known in the art, or after a 

sufficiently long period of time has passed which will erase the effects of spin labeling 
15 after the labeling procedure. After the effects of spin labeling are erased, the control 

magnetic field gradient and the control RF irradiation may be applied again. The control 

magnetic field gradient and control RF irradiation may be amplitude modulated 

according to a control modulation envelope. Furthermore, control RF irradiation may be 

frequency modulated according to the control modulation envelope. 
20 A second set of data may be acquired after the control procedure in step 410. As 

with the labeling dataset, it may later be converted to an image, if desired. Conversion of 

datasets to images may be performed in step 412. 

A perfusion image may be obtained in step 414 by subtracting the labeling image 

from the control image. Images or datasets (whichever appropriate) may be further 
25 processed by subtracting systematic residual errors in step 416. A particular kind of 

systematic error is described below in connection with Fig. 1 1 . 

One embodiment of the invention illustrated in Figs. 2-11 is not limited to the 

steps described above. Any of the steps may be repeated multiple times or in a different 

order than described. Furthermore, additional steps may be performed. A particular 
30 order and number of steps may be determined by one skilled in the art depending on the 

particular implementation and application of the invention. 
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Amplitude modulation envelopes for control and labeling are described in 
connection with Fig. 5. It should be noted that the scale displayed is displayed for 
illustrative purposes only and may not reflect the actual scale of modulation. 
Furthermore, the scale of modulation may vary for magnetic field gradient and for RF 
5 irradiation and from procedure to procedure. 

Control envelope 504 may be generally thought of as a square wave. However, 
in actual applications, the amplitude cannot change instantaneously, as is in the square 
wave, therefore envelope 504 is adapted to actual conditions. In general, any kind of 
envelope may be used, so long as a time scale for the amplitude modulation provided is 

10 fast as compared to the time T2 of the fluid and the time needed for the fluid to flow 
through the region where the effect of the labeling RF irradiation is large. Furthermore, 
any kind of wave may be used. In an alternative embodiment of the invention, envelope 
504 may resemble, for example, a triangular wave, or any other kind of wave, adapted to 
a particular procedure. For preferred results, control envelope 504 may have average 

15 amplitude equal to zero, in order to not perform any spin inversion. 

Control envelope 504 may be applied to both the magnetic field gradient and to 
control RF irradiation. Preferably, modulation may be synchronized between the 
magnetic field gradient and the RF irradiation. In an alternative embodiment of the 
invention, control amplitude modulation envelopes may differ for the magnetic field 

20 gradient and for the control RF irradiation, but may be similar to each other. The term 
"similar" may be used to refer to envelopes which are a slight variations of each other. 
In yet another embodiment of the invention, different modulation envelopes may be used 
for the magnetic field gradient and for the RF irradiation, so long as the physical 
properties of the net magnetization effects (Fig. 7) are generally similar to those in the 

25 first embodiment of the invention. In yet another embodiment of the invention, different 
magnetization results may be achieved. 

Labeling amplitude modulation envelope 502 may be an absolute value of control 
modulation envelope 504. In one embodiment of the invention, label amplitude 
modulation envelope 502 is an absolute value of a modified square function. An average 

30 labeling amplitude may be not zero in order to achieve spin labeling. In the embodiment 
of the invention described above, the average labeling magnitude may be positive, such 



- 12- 

as shown in Fig. 5. In an alternative embodiment of the invention, the average labeling 
magnitude may be negative. 

Labeling modulation envelope 502 need not be an absolute value of the modified 
square function. It may be any curve, as appropriate for a particular application. 
5 Furthermore, the labeling modulation envelope need not be an absolute value of the 
control modulation envelope - it may be a similar or an entirely different function, so 
long as the control RF irradiation mimics the magnetization transfer and other systemic 
effects of the labeling RF irradiation. Just as with the control modulation envelope, the 
labeling modulation envelope need not be identical for magnetic field gradient and the 
10 RF irradiation, although it is such in the embodiment of the invention illustrated in Figs. 
5 and 6. 

Fig. 6 is a diagram of a frequency modulation envelope for labeling and control 
RF irradiation. In one first embodiment of the invention previously described, frequency 
modulation envelopes 602 and 604 are similar to amplitude modulation envelops 502 

15 and 504, respectively. However, the invention is not limited to such envelopes. In an 
alternative embodiment of the invention, frequency modulation envelopes may be 
different from amplitude modulation envelopes. Furthermore, frequency modulation 
envelope 602 for labeling RF irradiation need not be an absolute value of the control RF 
irradiation envelope 604. In general, frequency and amplitude modulation envelopes 

20 502, 504, 602, and 504 may be varied depending on requirements of the particular 
application, as deemed appropriate by one skilled in the art. 

In one embodiment of the invention illustrated in Figs. 5 and 6, both amplitude 
and frequency modulation is fast compared to T2, in order to avoid decreases in 
efficiency due to transversal relaxation. In an alternative embodiment of the invention, 

25 modulation may be performed as fast as physically possible on the particular equipment 
used, or may be slower, depending on the fluids and tissues in question. 

Figure 7 illustrates the direction of magnetic fields with simultaneous RF and 
magnetic field gradient amplitude modulation and frequency modulation of the RF 
irradiation. 

30 When viewed in a reference frame spinning at the instantaneous frequency of the 

RF, a reference frame referred to as the rotating frame by those skilled in the art, the 
magnetization of the nuclear spins behaves as under the influence of an effective 
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magnetic field having two perpendicular vector components. The component along the 
direction of the field generated by the large, constant field applied within the MRI 
scanner is proportional to the difference between the frequency of the spins and the 
applied RF. The component along the perpendicular direction is equal to the amplitude 
5 of the applied RF. Magnetization tends to precess, or rotate, around the effective field. 
If the effective field direction changes only slowly, the precission of the magnetization 
will keep it approximately aligned with the effective field direction. The condition 
where the effective field direction changes only slowly is know to those skilled in the art 
as "adiabatic". 

10 The simultaneous amplitude and frequency modulation in accordance with 

embodiments of the invention will affect the RF amplitude and offset frequency in 
approximately equal proportion such that the direction of the effective field does not 
change. This preserves the adiabatic condition. When the amplitude modulation 
averages to a nonzero value, the magnetization will slowly follow the effective field 

15 direction. When the amplitude modulation averages to zero, precession in one direction 
will be canceled by precession in the opposite direction and there will be no net effect of 
the applied RF irradiation. 

Figs. 8 and 9 illustrate simulated results, while figures 10 and 1 1 illustrate results 
of the procedures performed on the human brain (in vivo). 

20 Simulations of inversion efficiency and direct off-resonance effects may be 

performed using the Bloch equation integration approach. Simulations may be 
performed for different phases of the modulations, for example, for 16 different phases 
of the modulation when the spin crosses the inversion plane. In simulation procedures, 
little variation of efficiency with phase was observed. Different velocities may be 

25 simulated and calculations of efficiency may be performed. 

Velocities from near zero to over lOOcm/s and higher have been simulated and 
labeling efficiency for a vessel was calculated assuming a parabolic flow profile. A 
square wave modulation was selected as the frequency and amplitude modulations 
envelope because the square wave modulation has the highest power efficiency for this 

30 particular application. An approximate square wave modulation at 1kHz has been used 
for the simulations. Real gradient systems cannot switch instantaneously, so the edges of 
the square wave transitions have been smoothed by appending half gaussians to the flat 
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portion of the modulation (Figs 5 and 6). In addition to the inversion efficiency, the 
direct effects of the applied fields to static tissue with a Tl time and a and T2 time 
typical of brain tissue have also been stimulated. 

Fig. 8 is a diagram of simulated residual subtraction error. Simulated residual 
5 subtraction error of B 1 of 1 8 mG and G of 0. 1 gauss/cm is shown as a function of 
distance from the labeling plane. For the parameters shown, subtraction error for the 
control approach is much less than the perfusion signal at a distance of 4cm from the 
labeling plane. 

Fig. 9 is a diagram of a simulated efficiency of one embodiment of the invention 

10 illustrated in Fig. 4. Simulated efficiency of greater than 90% has been achieved over a 
wide range of velocities even with RF amplitude reduced to 18mG (Fig. 9). The 
modulation approach of one embodiment of the invention is not affected by strong 
artifacts near the center frequency, has been shown to mimic the continuous flow driven 
adiabatic inversion efficiency with constant RF and gradient amplitude, in vivo. This 

15 embodiment removed systematic errors across a volume much larger than the constant 
RF and gradient approach. 

One embodiment of the invention has been implemented in- vivo on a General 
Electric VH/I 3 Tesla scanner. The invention is not limited to 3 Tesla scanner, and may 
be implemented on a wide range of MR apparatuses. Imaging has been performed with 

20 an echoplanar imaging sequence with a slice thickness of 7mm, a 24cm field of view, 

and an acquisition matrix of 64x64. The spin labeling was performed with the following 
timing: 1.3 second after a selective spatial saturation, the adiabatic RF irradiation was 
applied for 1.5 seconds and stopped 1.2 seconds before acquisition. Images were 
acquired with a TR of 5 seconds and each labeling and control image was averaged 12 

25 times for a total of 2 minutes of scanning. 

Fig. 10 is an illustration of in- vivo perfusion imaging achieved with one 
embodiment of the invention. Residual errors in the timing and shape of the RF and 
gradients caused a severe artifact on the subtraction images in certain specific locations 
(images 1 10). These and other errors may show up in other implementations of the 

30 invention. 

Such errors may be corrected by applying a small offset to the frequency of the 
RF irradiation. Such an offset may be constant and not affect the overall effectiveness of 
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the RF irradiation if the offset is relatively small compared to the overall frequency of 
the RF irradiation. In an alternative embodiment, the offset may be in a different 
proportion to the RF frequency. In yet another embodiment of the invention, a solution 
other than the frequency offset may be used to cure the undesired errors and unrelated 
5 effects. 

In an in- vivo implementation of one embodiment of the invention, a 100Hz offset 
to the frequency of the RF was used and it eliminated the artifacts of the undesired errors 
and produced images demonstrating strong perfusion contrast, as shown in images 1012. 
Additional errors may be present in the system, whether due to various 
10 imperfections, such as, for example, imperfections in the apparatus of a particular 

embodiment. Various mechanisms may be used to correct for some or all of those errors. 

Some residual errors may be detected if labeling and control are applied 
downstream to the region of interest being imaged. Fig. 1 1 is an illustration of 
comparative results of labeling upstream to an imaged region and downstream to the 
15 imaged region. 

For detection of residual errors, labeling of spins may be applied to a region 
downstream to the region of interest being imaged. In such a labeling, of course, no 
actual tissue perfusion is detected, because there may be no blood flow from the 
downstream region to the region of interest. Likewise, control RF irradiation may be 
20 applied to the downstream region. Theoretically, in case of downstream imaging, 

labeling and control datasets should be identical. However, in actual implementations 
that may be not the case due to residual and systematic errors. When labeling images are 
subtracted from control images, a small residual signal may still be present, as shown in 
images 1112. 

25 If desired, the data about residual or systematic errors may be subtracted from the 

perfusion image resulting from subtraction of labeling and control images. After these 
errors are removed, images illustrating blood perfusion may be received, as shown in 
images 1110. 

Additional image processing may be performed on the datasets received from the 
30 procedures. In the in- vivo procedures in the first embodiment of the invention, all 
recorded raw echo amplitudes were saved and transferred to a workstation for 
processing. Software written within the Instructive Data Language (IDL) environment 
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may be used to reconstruct the images. Various transformations may be applied to 
images, as appropriate for a particular implementation. Correction for image distortion 
and alternate k-space line errors may be performed on each image on the basis of data 
acquired during phase-encoded reference imaging. A side effect of this correction may 
5 be that image regions in which insufficient MR signal is present may be set to zero to 
allow for accurate measurement of the magnetic field. 

For labeling and control, magnitude images may be averaged. Gray and white 
matter may be automatically selected using the cerebrospinal fluid (CSF) suppression 
and CSF-white matter suppression. 

10 In the multi-slice study, CBF may be calculated from images acquired with and 

without arterial spin labeling. The averaged labeled images may be first subtracted from 
the averaged control images and then divided by the signal intensity on an image 
obtained in the absence of off-resonance saturation. 

The present invention is not limited to a particular implementation of the 

15 apparatus or dataset reconstruction system. It may be implemented on a variety of 

scanners and other MR equipment. Furthermore, it is not limited to MR imaging studies, 
and may be used for any MR study, as deemed appropriate by one skilled in the art. 

The invention may also be used for magnetic resonance angiography, the 
selective imaging of blood vessels by MRI. The invention may furthermore be used in 

20 measurement of flow in pipes or porous media and in a variety of other applications. 

The software implementation of the system is not limited to the implementation 
described above. Both system controls and an image reconstruction system may be 
implemented on a variety of computers or computer systems. 

A computer system for implementing an MR apparatus control and/or imaging 

25 system as a computer program typically includes a main unit connected to both an output 
device which displays information to a user and an input device which receives input 
from a user. The main unit generally includes a processor connected to a memory 
system via an interconnection mechanism. The input device and output device also are 
connected to the processor and memory system via the interconnection mechanism. 

30 It should be understood that one or more output devices may be connected to 

the computer system. Example output devices include a cathode ray tube (CRT) display, 
liquid crystal displays (LCD), printers, communication devices such as a modem, and 



- 17- 



audio output. It should also be understood that one or more input devices may be 
connected to the computer system. Example input devices include a keyboard, keypad, 
track ball, mouse, pen and tablet, communication device, and data input devices such as 
sensors. It should be understood the invention is not limited to the particular input or 
5 output devices used in combination with the computer system or to those described 
herein. 

The computer system may be a general purpose computer system which is 
programmable using a computer programming language, such as C++, Java, or other 
language, such as a scripting language or assembly language. The computer system may 

10 also include specially programmed, special purpose hardware. In a general purpose 

computer system, the processor is typically a commercially available processor, of which 
the series x86 and Pentium processors, available from Intel, and similar devices from 
AMD and Cyrix, the 680X0 series microprocessors available from Motorola, the 
PowerPC microprocessor from IBM and the Alpha-series processors from Digital 

15 Equipment Corporation, are examples. Many other processors are available. Such a 
microprocessor executes a program called an operating system, of which WindowsNT, 
UNIX, DOS, VMS and OS8 are examples, which controls the execution of other 
computer programs and provides scheduling, debugging, input/output control, 
accounting, compilation, storage assignment, data management and memory 

20 management, and communication control and related services. The processor and 
operating system define a computer platform for which application programs in 
high-level programming languages are written. 

A memory system typically includes a computer readable and writeable 
nonvolatile recording medium, of which a magnetic disk, a flash memory and tape are 

25 examples. The disk may be removable, known as a floppy disk, or permanent, known as 
a hard drive. A disk has a number of tracks in which signals are stored, typically in 
binary form, i.e., a form interpreted as a sequence of one and zeros. Such signals may 
define an application program to be executed by the microprocessor, or information 
stored on the disk to be processed by the application program. Typically, in operation, 

30 the processor causes data to be read from the nonvolatile recording medium into an 

integrated circuit memory element, which is typically a volatile, random access memory 
such as a dynamic random access memory (DRAM) or static memory (SRAM). The 
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integrated circuit memory element allows for faster access to the information by the 
processor than does the disk. The processor generally manipulates the data within the 
integrated circuit memory and then copies the data to the disk when processing is 
completed. A variety of mechanisms are known for managing data movement between 
5 the disk and the integrated circuit memory element, and the invention is not limited 
thereto. It should also be understood that the invention is not limited to a particular 
memory system. 

It should be understood the invention is not limited to a particular computer 
platform, particular processor, or particular high-level programming language. 

10 Additionally, the computer system may be a multiprocessor computer system or may 

include multiple computers connected over a computer network. It should be understood 
that each module (e.g. 100, 102a, 104a, 106, and 105) in FIG. 1A may be separate 
modules of a computer program, or may be separate computer programs. Such modules 
may be operable on separate computers. Data (e.g. 100 and 104a) may be stored in a 

15 memory system or transmitted between computer systems. The invention is not limited 
to any particular implementation using software or hardware or firmware, or any 
combination thereof. The various elements of the system, either individually or in 
combination, may be implemented as a computer program product tangibly embodied in 
a machine-readable storage device for execution by a computer processor. Various steps 

20 of the process may be performed by a computer processor executing a program tangibly 
embodied on a computer-readable medium to perform functions by operating on input 
and generating output. Computer programming languages suitable for implementing 
such a system include procedural programming languages, object-oriented programming 
languages, and combinations of the two. 

25 Having now described several embodiments, it should be apparent to those 

skilled in the art that the foregoing is merely illustrative and not limiting, having been 
presented by way of example only. Numerous modifications and other embodiments are 
within the scope of one of ordinary skill in the art and are contemplated as falling within 
the scope of the invention. 

30 What is claimed is: 



